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1.  SUMMARY 


We  have  investigated  the  kinetics  of  N2(X,  v)  quenching  by  silane.  These  studies 
were  a  first  step  in  understanding  the  kinetic  processes  important  in  silicon  nitride  coating 
systems  using  the  remote  plasma  enhanced  cttemical  vapor  deposition  technique.  Our 
experiments  were  carried  out  in  a  discharge-flow  reactor  and  used  metastable  energy  transfer 
diagnostics  to  detect  N2(X,  v)  in  vibrational  levels  from  zero  up  to  at  least  13. 

We  used  metastable  helium  Penning  ionization  to  detect  vibrational  levels  between 
zero  and  six,  and  N2(A)  energy  transfer  to  detect  levels  from  five  up  to  at  least  13. 
Metastable  helium  excites  N2(X,  v)  to  N2‘'‘(B^E‘'‘).  By  analyzing  the  resultant 
N2'''(B2e'''-X)  fluorescence  we  can  determine  tlte  vibrational  populations  of  N2(X,  v) 
quantitativelv.  N2(A^Eu)‘‘‘  excites  N2(X,  v)  to  N2(B^ng).  Analysis  of  the  subsequent 
N2(B^ng-A%u''')  fluorescence  gives  a  semiquantitative  estimate  of  the  vibrational 
distributions  for  levels  from  5  up  to  at  least  13. 


We  measured  the  rate  coefficients  by  monitoring  the  decrease  in  the  populations  of 
N2(X,  v)  as  a  function  of  [SiH4]  &nd  time  between  the  injection  of  the  silane  and  detection  of 
the  N2(X,  v).  We  analyzed  our  data  assuming  an  exponential  decay  of  the  N2(X,  v)  level 
under  investigation.  In  addition,  the  data  for  vibrational  levels  0  through  6  were  analyzed 
according  to  a  single-quantum  transfer  model.  The  exponential  decay  analysis  only  gives  the 
effective  removal  of  N2(X,  v)  in  the  observed  vibrational  level,  essentially  the  difference 
between  the  removal  of  vibrational  energy  from  a  given  level  to  those  below  it  and  the  feed 
into  the  vibrational  level  from  those  above  it.  The  single-quantum  transfer  model,  on  the 
other  hand  accounts  for  feed  from  higher  vibrational  levels.  Our  results  indicate  that  the  rate 
coefficient  for  quenching  N2(X,  v)  by  silane  varies  from  about  1  x  10'^^  cm^  molecule’^  s'* 
for  v=l  up  to  about  6  x  10^^  cm^  molecule'^  s'^  for  vibrational  levels  above  5. 


We  detected  fluorescence  from  SiH4(j'3)  at  4500  nm  when  silane  and  N2(X,  v)  were 
mixed,  thus  indicating  that  one  channel  for  the  quenching  reaction  is  vibration-vibration  (v-v) 
energy  transfer.  Comparing  the  relative  excitation  rate  coefficients  of  SiHj(j'3)  and  N20(j'3), 
excited  when  N2O  was  added  to  a  flow  of  N2(X,v),  indicated  that  only  about  15%  of  the 
total  N2(X,  v)  quenching  results  in  SiH4(i'3)  excitation.  This  assumes  that  N20(*'3) 
excitation  is  the  only  channel  for  N2(X,  v)  quenching  by  N2O.  Presumably  the  rest  of  the 
quenching  appears  in  other  vibrational  modes  of  SiH4,  although  >'3  excitation  would  be 
expected  to  be  favored  since  v-v  transfer  to  is  resonant  to  with^  kT  for  single  quantum 
transfer  from  all  N2(X,  v)  vibrational  levels  Itelow  16. 

We  had  planned  to  investigate  the  possibility  that  N2(X,  v)  quenching  would  lead  to 
SiH4  dissociation  by  monitoring  H-atom  production  using  vacuum  ultraviolet  (vuv)  resonance 
fluorescence  at  121.6  nm.  Unfortunately,  our  vuv  detection  sensitivity  below  130  nm  was  so 
small  H-atom  resonance  fluorescence  was  undetectable  even  when  significant  number 
densities  of  H  atoms  were  purposely  added  to  the  reactor.  This  degeneration  seems  to  be  a 
recent  problem  since  our  system  sensitivity  was  quite  good  several  years  ago. 
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We  began  a  series  of  experiments  in  which  we  coated  single  crystal  silicon  chips  with 
silicon  nitride.  These  experiments  were  performed  in  another  discharge-flow  reactor.  In  this 
apparatus,  active  nitrogen  was  generated  in  a  microwave  discharge  sustained  in  a  flowing 
mixture  of  nitrogen  in  argon.  Silane  was  added  to  the  discharge  afterglow  and  the  reaction 
products  impinged  upon  a  heated  substrate  further  downstream.  Infrared  absorption 
spectroscopy  showed  that  the  coatings  generated  were  indeed  silicon  nitride,  although  some 
hydrogen  was  incorporated  into  the  films.  In  addition  to  the  coatings,  the  active  nitrogen/ 
silane  interaction  generated  a  quantity  of  fuie  (micron  sized  particles),  white  powder. 

Infrared  absorption  spectroscopy  indicated  this  powder  to  be  silicon  nitride. 
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2.  INTRODUCTION 


Silicon  nitride  films  find  wide  application  as  abrasion-,  corrosion-,  and  oxidation- 
resistant  coatings  to  protect  machined  parts  in  high  temperature  environments.  They  are  used 
in  numerous  ways  including  radiant  heat  shields,  rocket  nozzle  inserts  and  crucibles  for 
melting  and  refuting.  Some  applications  require  relatively  cool  substrate  temperatures  during 
the  coating  process  so  that  other  species  on  the  substrate  don’t  evaporate  or  migrate  to 
undesirable  areas  of  the  substrate.  The  limited  mechanical  strength  of  some  substrate 
structures  necessitates  cool  substrate  temperatures  during  the  coating  process.  ^ 

Switching  from  normal  chemical  vapor  deposition  (CVD),  for  which  substrate 
temperatures  typically  are  700  to  1000**C,  to  plasma  enhanced  chemical  vapor  deposition 
(PECVD)  allows  coating  down  to  temperatures  on  the  order  of  200  to  500“C.^  Typical  con¬ 
ditions  for  PECVD  involve  placing  the  substrate  to  be  coated  inside  a  discharge  maintained 
in  mixture  of  SiH4  and  either  N2  or  NH3.  This  procedure  typically  results  in  the 
incorporation  of  between  10  and  30  atom  percent  hydrogen  in  the  coating.  For  some 
applications,  this  level  of  hydrogen  is  unacceptable. 

A  variation  of  PECVD,  which  also  is  suitable  for  substrate  temperatures  on  the  order 
of  200  to  SOO^’C,  and  which  results  in  undetectable  concentrations  of  Si-H  or  N-H  bonds  in 
deposited  films,  has  been  dubbed  by  one  of  its  developers  "remote  plasma  enhanced  chemical 
vapor  deposition"  (RPECVD).^  In  this  technique,  the  discharge  is  removed  from  the 
substrate.  Silane  is  added  to  the  effluents  of  a  discharge  in  flowing  nitrogen.  Energetic 
species  in  the  discharged  nitrogen  (generally  called  "active  nitrogen")  decompose  the  silane; 
the  reaction  products  from  this  decomposition  then  impinge  upqn  the  substrate  and  lay  down 
a  film  of  silicon  nitride. 

While  both  PECVD  and  RPECVD  have  been  studied  and  used  for  a  number  of  years, 
the  important  reactions  occurring  in  the  plasma  and  plasma  afterglow  are  not  well  under¬ 
stood.  This  lack  of  understanding  hampers  efforts  to  improve  the  efficiency  and  purity  of 
silicon  nitride  coating  processes. 

Several  groups  have  suggested  that  atomic  nitrogen  is  of  prime  importance  in  the  dis¬ 
sociation  of  silane.^’^  Indeed,  one  study  has  shown  that  the  extent  of  silane  decomposition 
in  the  active  nitrogen  scaled  linearly  with  the  number  density  of  nitrogen  atoms  in  the  active 
nitrogen  flow.^  Atomic  nitrogen  cannot  be  the  only  species  controlling  silane  decomposition 
because  the  reaction  between  nitrogen  atoms  and  silane  is  endoergic  by  about  half  of  an 
electron  volt.^  Furthermore  we  have  shown^  that  the  rate  coefficient  for  the  reaction 
between  atomic  nitrogen  and  silane  is  insignificant  (k  <  8  x  10'^'*  cm^  molecule'^  s'*).  The 
reactions  of  nitrogen  atoms  with  SiH^  radicals  with  x  ^  3  are  exoergic  and  undoubtedly  do 
contribute  to  the  overall  decomposition  of  silane.  Some  species  other  than  atomic  nitrogen, 
however,  must  be  responsible  for  the  initial  decomposition  step. 
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While  some  investigators  think  that  vibrationally  excited  nitrogen  is  crucially 
important  in  the  initial  decomposition  step^'^,  several  other  research  groups  have  concluded 
that  vibrationally  excited  nitrogen  does  not  appear  to  be  important  in  the  silane 
decomposition  process.^’^’^ 

Vibrationally  excited  nitrogen  could  contribute  to  silane  dissociation  in  two  different 
ways.  In  both  paths  the  excited  nitrogen  would  transfer  several  quanta  of  vibrational  energy 
to  the  vi  or  modes  of  silaite, 

N2(v)  >  SiH4  -  N2(v-1)  >  SiH4(»'3,»'j)  .  (D 


These  two  modes  of  silane  are  essentially  resonant  (within  one  kT)  for  one  quantum  transfers 
from  nitrogen  molecules  that  are  excited  as  high  as  the  fifteenth  vibrational  level  ( -  4  eV). 
Since  resonant  energy  transfer  procesMS  generally  proceed  facilely  this  vibrational  energy 
transfer  process  should  be  efficient. 

The  transfer  of  two  quanta  from  N2(v)  to  SiH4  is  sufficient  to  render  the  reaction 
between  ground-state  nitrogen  atoms  and  silane  exoergic, 

SiH4(2  1^3,  vj)  N('^S)  -  SiH3  +  NH  +  0.12  eV  .  (2) 

The  SiH^  fragments  can  then  be  decomposed  further  in  a  series  of  N-atom  reactions  similar 
to  reaction  2. 

The  other  possibility  is  for  the  nitrogen  to  transfer  vibrational  energy  to  the  silane 
until  the  silane  obtains  enough  internal  energy  to  dissociate  unimolecularly, 

SiH4(n  1^3, I'll  -*  SiH2  +  H2  (n>9)  ,  (3a) 

or 

SiH4(n  SiH3  +  H  (n^l5)  .  (3b) 


The  SiH4  decomposition  would  require  at  least  nine  quanta  of  nitrogen  vibrational  energy  to 
effect  the  lowest  energy  decomposition  step,  that  into  SiH2  and  H2.  If  reaction  3b  is  the 
preferred  decomposition  pathway,  15  N2(X,  v)  quanta  will  be  required.  The  two  branches 
can  be  distinguished  rather  easily  because  reaction  3b  will  produce  atomic  hydrogen,  whereas 
reaction  3a  will  not.  Vacuum  ultraviolet  resonance  fluorescence  measurements  provide  a 
simple  yet  sensitive  diagnostic  for  hydrogen  atoms. 

The  vibrational  pumping  of  the  silane  could  proceed  directly  from  a  series  of  energy 
transfer  collisions  between  the  vibrationally  excited  nitrogen  and  the  silane. 
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N2(v)  +  SiH4(m 


N2(v-1)  +  SiH4(m+l  vyv{)  ; 


(4) 


alternatively  a  number  of  energy  pooling  reactions  between  vibrationally  excited  silane 
molecules  could  occur: 

» 

SiH4(n  +  SiH4(m  -*  SiH4(n-l  >'3,>'i)  +  SiH4(m+l  V2,vi)  ,  (5) 


where  m>n. 

O’Keefe  and  Lampe  invoked  energy-pooling  processes  to  explain  their  observations  of 
silane  decomposition  in  mixtures  of  SiF4  and  SIH4  that  had  been  irradiated  by  a  CO2  laser 
at  wavelengths  that  are  absorbed  only  by  the  SiF4.  They  inferred  that  the  laser  irradiation 
resulted  in  vibrational  excitation  of  the  SiF4,  which  subsequently  transferred  its  energy  to  the 
SiH4.  The  vibrationally  excited  SiH4  then  participated  in  a  number  of  energy  pooling 
collisions  that  ultimately  led  to  decomposition  into  SiH2  and  H2.  There  was  no 
decomposition  of  the  SiF4,  only  the  SiH4. 

The  temporal  behavior  of  the  decomposition  following  a  laser  pulse  was  quite  inter¬ 
esting.  Although  the  laser  pulse  lasted  only  about  2  /is,  silane  decomposition  continued  for 
roughly  a  millisecond.  The  decomposition  rate  reached  a  peak  about  200  ns  after  the  laser 
pulse  and  then  decayed  with  a  lifetime  of  approximately  300  The  most  intriguing  feature 
is  that  the  decomposition  process  appeared  to  have  an  incubation  period.  The  rate  of 
decomposition  was  an  order  of  magnitude  greater  during  the  second  100  /xs  than  during  the 
first  100  ^s  following  the  laser  pulse. 

The  incubation  period  prior  to  rapid  silane  decomposition  that  O’Keefe  and  Lamp 
noted  is  reminiscent  of  observations  in  two  recent  flowing  afterglow  studies  on  the  decompo¬ 
sition  of  silane  in  active  nitrogen.^'^  Both  studies  noted  an  incubation  period  prior  to  silane 
decomposition.  Horie  et  al.^  observed  that  both  silane  and  atomic  nitrogen  number  densities 
remained  relatively  constant  over  periods  of  time  up  to  35  ms,  depending  upon  the  ratio  of 
initial  number  densities  of  atomic  nitrogen  to  silane,  and  then  decreased  fairly  rapidly  after 
this  incubation  period.  For  example,  under  one  set  of  conditions  that  they  investigated,  the 
silane  and  N-atom  number  densities  decreased  very  little  for  more  than  20  ms  after  the  silaiK 
and  active  nitrogen  had  been  mixed.  During  the  next  6  to  8  ms,  however,  the  silane  and  the 
N  atoms  were  consumed  completely.  This  peculiar  behavior  could  be  explained  if  something 
like  20  ms  were  required  for  nitrogen  to  couple  adequate  internal  energy  to  the  silane  before 
it  could  decompose.  Then  once  the  decomposition  was  begun,  its  rate  could  accelerate 
rapidly  if  some  reaction  chain  continued  further  decomposition. 

DeJoseph  proposed  a  mechanism  that  gave  an  alternative  explanation  to  the  incubation 
time.^  He  posmlated  that  the  primary  reaction  necessary  to  begin  silane  decomposition  was, 

N2(X,  v^lO)  +  SiH4  -*  N2  +  SiH2  +  H2  .  (6) 
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The  SiH2  was  then' further  decomposed  by  reactions  with  nitrogen  atoms.  The  novel  idea  in 
this  postulated  decomposition  mechanism  was  a  series  of  three  reactions  that  regenerated  the 
N2(v): 


SiN  +  N  -  Si  +  N2(X,  v^20)  , 

(7) 

Si  +  N2  +  M  -  SiN2  +  M  , 

(8) 

SiN2  +  N  -  SiN  N2(X,  v^13) 

(9) 

The  net  effect  of  these  three  reactions  is  to  generate  two  N2(v)  molecules  in  the  recombi¬ 
nation  of  two  N' atoms.  As  the  reaction  progresses,  therefore,  the  number  density  of  the 
N2(X,  V  ^  10)  builds  up  and  this  build  up  accelerates  the  rate  of  silane  decomposition.  With 
this  set  of  reactions  at  its  core,  DeJoseph  constructed  a  kinetic  model  for  silane 
decomposition  that  reproduced  moderately  well  the  temporal  observations  of  Horie  et  al. 
using  a  reasonable  set  of  rate  coefficients. 

The  basic  mechanism  for  silane  decomposition  described  by  reactions  6  to  9  can  be 
distinguished  experimentally  from  that  which  incorporates  reactions  1  to  5  because  the 
temporal  variation  of  the  number  densities  of  N'2(X,  v  ^  10)  will  differ  dramatically  between 
the  two  models.  The  first  mechanism,  reaction  1  to  5,  will  have  the  N2(X,  v^  10)  number 
densities  diminish  over  time.  Expectations  based  upon  the  second  mechanism,  reaction  6  to 
9,  on  the  other  hand,  are  that  the  N2(X,  v  ^  10)  number  densities  will  increase  by  several 
orders  of  magnitude  over  time.  Identifying  the  correct  mechanism  requires  a  diagnostic  for 
vibrational  levels  of  ground-state  nitrogen  greater  than  ten.  We  developed  this  capability 
recently  during  a  course  of  investigations  on  the  kinetics  of  N2(X,  v)  in  prospective  chemical 
laser  systems. 

The  basic  goals  of  this  program  were  to  try  to  determine  the  importance  vibrationally 
excited  nitrogen  plays  in  silane  decomposition  in  silicon  nitride  coating  systems.  Our 
primary  focus  was  on  characterizing  tte  reactions  between  N2(X,  v)  and  SiH4.  To  this  end 
we  studied  the  kinetics  of  N2(X,  v)  quenching  by  SiH4  in  some  detail.  In  additon  to 
determining  rate  coefficients  as  a  function  of  vibrational  level,  we  showed  that  the  quenching 
process  can  be  characterized  as  a  series  of  single-quantum  transfers  and  that  one  channel  for 
the  reaction  is  SiH4(i'3)  excitation. 

We  also  began  a  series  of  experiments  aimed  at  forming  silicon  nitride  coatings  on 
single-crystal  substrates.  Although  these  experiments  were  limited  in  extent,  successful 
coatings  were  generated  with  little  difficulty.  In  addition  silicon  nitride  powder  was  also 
formed  in  these  experiments. 
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3.  EXPERIMENTAL 


The  experiments  were  all  done  in  a  discharge-flow  reactor.  Different  configurations 
were  used  depending  upon  the  effects  being  studied.  The  bulk  of  the  experiments,  designed 
to  measure  rate  coefficients  for  N2(X,  v)  removal  by  SiH4,  used  the  apparatus  shown 
schematically  in  Figure  1.  It  consists  of  a  4.6-cm  i.d.  quartz  tube,  50  cm  long.  N2(X,  v)  is 
injected  into  the  upstream  end  of  the  reactor,  and  is  detected  at  the  downstream  end  by 
monitoring  fluorescence  excited  by  the  interaction  of  N2(X,  v)  with  either  N2(A)  or 
He(2  ^S).  A  loop  injector  that  can  be  moved  along  the  length  of  the  reactor  mixes  silane 
with  N2(X,  v).  The  flow  reactor  was  pumped  by  a  roots  blower/forepump  combination 
capable  of  producing  flow  velocities  of  5  x  10^  cm  s'*  at  a  pressure  of  1  Torr.  Generally 
much  more  modest  flow  velocities  were  required  and  the  flow  reactor  access  to  the  pump 
was  throttled  by  a  downstream  butterfly  valve. 


Figure  t.  Discharge  flow  reactor  for  studying  the  kinetics  of  reactions  involving  N2(X,v) 

N2(X,  v)  is  prepared  in  a  microwave  discharge  through  a  mixture  of  nitrogen  dilute  in 
helium.  The  effluents  of  the  discharged  mixture  then  pass  through  a  nickel  screen  prior  to 
their  entry  into  the  flow  reactor.  The  screen  removes  the  atoms  and  deactivates 
electronically  excited  metastables,  but  has  relatively  little  effect  on  vibrationally  excited  N2. 
Our  previous  observations*^  indicate  that  the  Ni  screen  reduces  N  atom  number  densities 
about  a  factor  of  SO,  to  a  level  of  about  1  x  10**  cm^  molecule'*  s'*.  This  reduction  is,  of 
course,  dependent  to  some  extent  upon  the  screen  mesh  and  contact  time  between  the 
discharge  effluents  and  the  screen.  In  our  experiments  the  screen  is  about  24  mesh  and 
consists  of  two  coaxial  cylinders  about  5  cm  long  inserted  into  the  12.7-mm  o.d.  discharge 
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tube  downstream  from  the  discharge.  They  cover  the  opening  in  the  side  of  the  discharge 
tube  through  which  the  effluents  must  pass  to  enter  the  flow  reactor. 

Metastable  He(2  ^S)  and  N2(A)  enter  tl^  flow  reactor  through  a  10  cm  o.d.  hook¬ 
shaped  injector  inserted  into  the  flow  reactor  at  its  downstream  end.  The  metastables  are 
generated  in  the  upstream  end  of  the  injector  and  emanate  from  the  injector  along  the  flow 
reactor  axis.  A  hollow  cathode,  dc  discharge,  biased  at  290V  dc,  generates  He(2  ^S)  if  the 
injector  flow  is  pure  helium,  or  metastable  Ar(^PQ  2)  if  a  few  percent  argon  is  added  to  the 
helium  flow.  The  metastable  argon  is  used  to  produce  the  N2(A)  metastables  by  excitation 
transfer.  In  this  case,  nitrogen  is  added  in  die  injector  just  downstream  from  the  discharge 
cathode.  The  lack  of  emission  from  ^^(C)  or  N^fB,  v^2)  when  the  He  metastables  are 
added  to  cold  N2  shows  that  both  He'*'  and  He^  are  absent  from  the  flow  of  metastables 
entering  the  flow  reactor. 

Gas  flow  rates  are  measured  by  electronic  mass  flow  meters,  or  in  a  few  instances  by 
rotameters,  that  have  been  calibrated  by  measuring  the  rates  of  pressure  increase  as  a 
function  of  time  when  gases  flow  into  a  calibrated  volume.  Because  some  of  our  initial 
results  were  puzzling,  the  mass  flow  meters,  were  extensively  recalibrated.  The  only 
significant  change  from  this  procedure  was  that  the  SiH4  calibration  was  about  15%  less  than 
we  had  anticipated  based  upon  the  flow  meter  calibration  with  argon  and  the  standard 
correction  procedure  involving  the  ratio  of  gas  specific  heats.  Capacitance  manometers 
measured  the  flow  reactor  pressures. 

A  0.5-m  monochromator,  sensitive  to  radiation  between  200  and  900  ran,  detected 
radiation  in  the  flow  reactor;  the  N2  first-negative  bands,  N2(B  -  X  between  350 
and  500  ran,  from  the  He*(2  ^S)/N2(X,  v^6)  interaction;  the  N2  first-positive  bands, 

N2(B  ^Ilg  -  A  between  560  and  900  ran,  from  the  N2(A  ^Eu)/N2(X,  v^5)  interaction; 
and  the  N2  Vegard-Kaplan  bands,  N2(A  ^Ey  -  X  *Ep  between  250  and  400  nm,  diagnostic 
of  N2(A)  number  densities.  The  monochromator  had  a  1200  groove  mm‘*  grating  blazed  at 
250  nm  and  a  thermoelectrically  cooled  photomultiplier.  Signals  from  the  photomultiplier 
were  processed  by  either  a  photon-counting  rate  meter  or  a  picoammeter.  Their  outputs  were 
recorded  by  a  strip-chart  recorder  and  stored  digitally  by  a  lab  computer  system.  A  least- 
squares,  spectral  fitting  procedure,  described  in  some  detail  previously,  determined  band 
intensities  from  the  spectra. 

Another  set  of  experiments  used  the  apparatus  shown  schematically  in  Figure  2. 

These  experiments  were  designed  to  monitor  infrared  fluorescence  from  SiH4  excited  by 
vibration-vibration  exchange  from  N2(X,  v).  For  these  experiments  a  5-cm  i.d.  flow  mbe 
was  configured  to  allow  spectral  observations  along  the  axis  of  the  flow  reactor.  The 
upstream  end  of  the  flow  reactor  was  sealed  with  a  CaF2  window  and  placed  directly  in  front 
of  the  monochromator  slit.  The  effective  field  of  view  was  about  32  cm  long.  Generally  the 
N2(X,  v)  entered  directly  in  front  of  the  window,  normal  to  the  axis  of  the  field  of  view. 
Silane,  or  in  some  cases  nitrous  oxide,  was  added  through  a  loop  injector  situated  directly  in 
front  of  the  window.  In  some  experiments,  the  N2(X,  v)  and  silane  flowed  through  a  50-cm 
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Figure  2.  Schematic  of  apparatus  for  observing  SiH4(('3)  infrared  emission 


long,  2.5-cm  i.d.  side  arm  before  entering  tlK  S-cm  tube  within  the  monochromator  field  of 
view. 


Most  infrared  spectral  observations  were  made  between  2000  and  5000  nm  using  the 
O.S-m  monochromator  equipped  with  a  150  groove  mm'^  grating  blazed  at  3000  nm  and  a 
liquid  nitrogen  cooled  InSb  detector.  In  a  few  cases  a  300  groove  mm'^  grating  blazed  at 
1000  nm  and  a  liquid  nitrogen  cooled  intrinsic  Ge  detector  were  used  with  the 
monochromator  to  allow  the  spectral  region  between  1000  and  1600  nm  to  be  monitored.  In 
both  cases,  light  from  the  observation  region  was  chopped  prior  to  entering  the 
monochromator  and  signals  from  the  detectors  were  processed  with  a  phase-sensitive 
amplifier.  The  ouqnit  from  the  lock-in  amplifrer  was  recorded  by  a  strip-chart  recorder  and, 
in  addition,  stored  digitally  in  a  computer  file. 

A  third  set  of  experiments  were  directed  at  coating  several  substrates  with  silicon 
nitride  films.  These  experiments  used  a  different  5  cm  diameter,  quartz  flow  reactor, 
illustrated  schematically  in  Figure  3.  It  was  pumped  by  a  rmary  mechanical  pump.  Active 
nitrogen  was  generated  at  the  upstream  end  of  the  flow  reactor  in  a  microwave  discharge  of 
nitrogen  dilute  in  argon.  Silane  was  injected  through  a  loop  injector  whose  position  could  be 
varied  along  the  length  of  the  reactor,  and  the  aaive  nitrogen,  silane  mixture  impinged  upon 
a  substrate  mounted  on  a  copper  cap  on  the  end  of  a  rod  heater.  Some  initial  experiments 
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Figure  3.  Schematic  of  apparatus  for  coating  substrates  with  silicon  nitride 


used  sapphire  substrates,  but  later  experiments  used  fragments  of  single  crystal  silicon 
wafers,  polished  on  both  sides.  The  substrate  was  sandwiched  between  the  cap  on  the  rod 
heater  and  a  copper  ring  that  was  screwed  into  the  cap.  In  addition  a  thermocouple  was 
wedged  under  (tic  copper  retainer.  The  output  of  the  thermocouple  was  wired  to  a 
commercial  temperature  controller.  As  in  all  other  experiments,  flow  rates  were  monitored 
by  electronic  mass  flow  meters.  We  found  it  necessary  to  use  a  mass  flow  controller  to 
maintain  a  constant  flow  of  silane.  The  other  gas  flows  remained  constant  so  long  as  their 
gas  cylinder  delivery  pressures  and  needle  valve  settings  remained  constant.  The  pressure  in 
this  reactor  was  measured  with  a  capacitance  manometer.  Typical  conditions  included  flow 
rates  of  argon  carrier  gas  and  nitrogen  of  880  and  90  micromoles  s'^  respectively,  and  a 
total  pressure  of  3  Torr. 
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4.  ENERGY-TRANSFER  BASED  DIAGNOSTICS  OF  NjCX,  t) 


We  have  developed  techniques  recently  for  determining  number  densities  of  vibra- 
tionally  excited  nitrogen  in  vibrational  levels  up  through  about  fifteen.  Two 

different  methods  are  used  to  probe  different  ranges  of  vibrational  excitation.  The  first, 
which  gives  absolute  N2(X,  v)  number  densities  for  vibrational  levels  up  through  about  six, 
involves  measuring  the  vibrational  distribution  of  the  nitrogen  first-negative  bands  when  they 
are  excited  in  the  Penning-ionization  reaction  between  metastable  helium  atoms  and 

N2(X,  V), 


He  *(2^S)+  N2(X,v)  -  He(*S)  +  N2(B,v).  (10) 

The  ground-electronic-state  vibrational  distribution  of  nitrogen  can  be  deduced  from  the 
N2(B)  vibrational  distribution  and  the  Franck-Condon  factors  that  connect  the  ground 
electronic  state  of  the  neutral  molecule  to  the  electronically  excited  state  of  the  molecular  ion 
(vide  infra). 

The  second  N2(X,  v)  diagnostic  gives  semi-quantitative  estimates  of  N2(X,  v)  number 
densities  for  vibrational  levels  greater  than  4  (about  S  to  Metastable  nitrogen 

molecules,  N2(A^Eu),  excite  N2(X,  v-5-15)  to  N2(B,  vi:l2)  quite  efficiently, 


N2(A^  Eu)  +  N2(X,v)  N2(X)  *  N2(B\,v^12),  dD 

where  k2  4  x  10'^^  cm^  molecules'^  s'^.  Measurements  of  the  intensity  and  vibrational 
distribution  of  the  N2(B)  fluorescence  in  the  presence  of  known  amounts  of  N2(A)  will  give 
the  number  densities  of  these  highly  excited  levels  of  N2(X,  v)  to  within  factors  of  two  to 
three.  Relative  variations  in  N2(X,  v)  number  densities  can  be  determined  to  within  about 
20%. 


We  have  employed  the  second  N2(X,  v)  diagnostic  to  determine  rate  coefficients  for 
vibrational  energy  transfer  from  N2(X,  v)  to  a  number  of  other  molecules.  The  technique 
involves  comparing  the  N2(B)  intensities  from  N2(A)  excitation  of  N2(X,  v)  between 
conditions  un^r  which  the  mixing  time  of  the  quencher  and  the  vibrationally  excited 
nitrogen  is  quite  short,  2  to  3  ms,  and  conditions  under  which  the  contact  time  is  much 
longer,  20  to  40  ms.  This  approach  allows  one  to  distinguish  between  reductions  in  N2(B) 
intensity  caused  by  electronic  quenching  of  N2(B)  from  those  resulting  from  a  reduction  in 
the  N2(X,  v)  number  density. 

These  measurements  use  a  relatively  clean  source  of  N2(X,  v),  one  that  is  essentially 
free  from  accompanying  atomic  nitrogen  or  electronically  excited  metastables.  This  source 
involves  passing  the  effluents  from  a  discharge,  in  a  mixture  of  nitrogen  in  helium,  through  a 
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nickel  screen.  The  nickel  screen  recbices  atomic  nitrogen  number  densities  in  the  afterglow 
by  factors  of  SO  or  more,  but  attenuates  the  N2(X,  v)  very  little. 

4.1  Praning  lonizatioa 

Our  Penning-Ionization  diagnostic  for  vibrationally  excited  nitrogen  in  discharge 
afterglows  follows  along  the  lines  of  the  pioneering  work  of  Schmeltekt^f  et  al.^^  and  Young 
and  Hom.^’^^  Mixing  metastable  helium  atoms  with  a  flow  of  molecule  nitrogen  results  in 
strong  emission  of  the  nitrogen  first-negative  system,  Since  the 

Penning-Ionization  process  follows  a  Franck-Condon  excitation  pathway,  the  vibrational 
distribution  in  the  neutral,  ground  state  will  determine  the  distribution  observed  in  the  upper, 
ionic  state.  One  problem  with  this  approach  is  that  care  must  be  taken  n(X  to  have  any  He*^ 
or  He2  in  the  flow  of  metastable  helium.  Both  of  those  species  also  excite  N2(B)  quite 
strongly  in  charge-transfer  reactions,  but  with  an  N  2(B)  vibrational  distribution  that  is 
decidedly  non-Franck-Condon.  We  have  been  careful  to  eliminate  ions  in  our  system. 

The  rate  of  Penning-ionization  between  metastable  helium  atoms  and  molecular 
nitrogen  is  proportional  to  the  Franck-Condon  factors  that  connect  the  N2(B)  and  N2(X) 
states.  One  can  calculate  the  vibrational  distrilxition  in  the  final  state,  therefore,  knowing 
only  the  vibrational  distribution  in  the  lower  state  and  the  Franck-Condon  factors  that  couple 
the  two  states.  Thus, 


N, 


V' 


Ny^Qv’^ 


(12) 


where  v'  and  v"  represent  the  vibrational  levels  of  the  upper  and  lower  states,  respectively, 
and  qy'yo  is  the  Franck-Condon  factor  coupling  them.  We  have  calculated  a  comprehensive 
set  of  Franck-Condon  factors  for  N  2(B)  -  N2(X)  transitions,  and  have  tabulated  them 
elsewhere.^^ 

We  determine  relative  N2(X,  v)  distributions  by  fmding  the  model  distribution  that 
can  be  used  in  Eq.  (12)  to  predict  N2(B)  vibrational  distributions  that  are  most  like  those 
observed  euerimentally.  The  N2(X,  v)  vibrational  distributions  usually  follow  a  modified 
Treanot^*^'  model,  which  simplifies  the  fitting  process.  We  do  have  fitting  procedures, 
however,  to  determine  populations  for  cases  for  which  the  Tteanor  model  is 
inapplicable. 

The  modified  Treanor  distribution  has  been  outlined  by  Caledonia  and  Center,^^  and 
Dilonardo  and  Capitelli.^  For  low  vibrational  levels  the  distribution  is  that  given  by 
Treanor  et  al.:^ 
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where  6i  is  the  Boltzmann  vibrational  ten^rature  referenced  to  v”  =  1,  T  is  the  ambient 
gas  translational  temperature,  and  and  are  ^pectroscq>ic  constants  in  units  of  cm*^ 
The  Boltzmann  vibrational  temperature  is  given  by 
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k  ln(Ny»al/Ny»3Q) 


(14) 


where  k  is  Boltzmann’s  constant. 

This  distribution  goes  through  a  minimum,  called  the  Treanor  minimum,  at  v*.  given 
by 
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For  vibrational  levels  above  the  Treanor  minimum,  the  Product  vNy  is  essentially  constant.^4 
The  resulting  distribution  for  v  >  v*  is 

Ny.  _(v  •  -1)  exp  (-l.4388(v 

Ny-.0  V 

where  the  various  parameters  in  Eq.  (16)  are  determined  by  requiring  the  two  distributions, 
i.e.,  Eqs.  (13)  and  (16),  be  equal  at  v"  and  v*-l.  Rather  than  use  an  explicit  Gj  value,  we 
base  our  distributions  on  an  effective  Gj  which  best  Hts  all  the  data  points. 

Fitting  the  Penning-ionization  spectra  to  a  ground-state  vibrational  model,  as  just 
described,  gives  only  the  relative  distribution  among  ground-state  vibrational  levels.  We 
place  this  relative  distribution  on  an  absolute  basis,  however,  by  determining  the  absolute 
number  density  of  N2(X,  v”  =0).  This  is  accomplished  by  noting  changes  in  intensity  of 
N2  (B,  v'=  0,1)  when  the  N2(X,  v)  producing  discharge  is  struck  relative  to  the  intensities 
when  the  discharge  is  off.  We  know,  of  course,  the  absolute  N2(X,v'’=0)  number  density 
when  the  discharge  is  off. 

4.2  N2(A)  Energy  Transfo* 

The  Penning-ionization  diagnostic  is  inadequate  for  studies  on  N2(X,v'’^6)  as 
discussed  earlier.  To  overcome  this  problem  we  developed  a  semiquantitative  technique  for 
determining  number  densities  of  N2(X,  v”  25).  Our  diagnostic  method  for  monitoring  N2(X, 
v"  25)  is  based  upon  measuring  nitrogen  flrst-positive  emission  intensities  excited  in  the 
reaction  between  N2(A)  and  N2(X,  v").  These  intensities  are  directly  proportional  to  N2(X, 
v")  number  densities.  The  observed  intensities  must  be  corrected  for  any  quenchii^  of  the 
first-positive  emission  before  they  can  be  related  to  N2(X,  v")  number  densities. 
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The  iittensity  of  the  N2  first-positive  emisssion  is  given  by 


I  =  k,[N2(B)]  = 


kii[N2(X.v'’)l[N2(A)l 


kQi[Qil 

kr 


(17) 


where  kgi  is  the  rate  coefficient  for  N2(B)  quenching  by  quenchers  Qj  and  k,  is  the  N2(B) 
radiative  decay  rate.^^  The  buffer  gas  {wessure.  primarily  helium  is  ^  1  torr.  This 
minimires  quenching  by  the  badi  gas  species  in  the  flow  tube.^^  Rearranging  Eq.  (17)  gives 


[N2(X,v'’)] 


k,(N2(B)l  ,  ^  y  koiEQi] 
knlN2(A)l  ^  i  k. 


(18) 


While  determining  absolute  values  for  [N2(A)]  or  [N2(B)]  is  difficult  and  therefore 
quite  uncertain,  measuring  the  ratio  of  these  two  number  densities  accurately  is  straight¬ 
forward  because  only  relative  sensitivities  of  the  detection  system  at  the  various  wavelengths 
are  required. 


Although  primarily  a  global  monitor  of  N2(X,  v"),  this  diagnostic  does  provide  some 
state  specificiQr.  The  hi^r  vibronic  levels  of  N2(B)  are  accessible  only  from  N2(X,  v") 
vibrational  levels  at  the  high  end  of  the  diagnostic  range.  For  example  N2(B,v'  =  10)  can  be 
excited  only  from  N2(X,  v"  2;  12).  Thus,  we  track  the  varying  behavior  of  N2(X,  v" 12) 
by  measuring  N2(B,v=10)  intensities. 
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5.  NjCX,  ▼)  QUENCHING  BY  SiH4 
5.1  ExperimeDtal  ApprosM^i 

Because  our  diagnostic  techniques  involve  adding  metastables  to  the  flow  reactor  and 
observing  the  fluorescence  excited  when  the  metastable  intoacts  with  SiH4,  we  need  to 
separate  diminutions  in  the  observed  fluorescence  caused  by  fluorescence  or  metastable 
(pienching  from  the  diminutions  reailting  frtnn  removal  of  N2(X,  v)  flcnn  die  reactor.  We 
separate  these  effects  by  varying  the  mixing  time  between  the  N2(X,  v)  and  the  SiH4  over  a 
wide  range. 

When  quencher  is  injected  immediately  behind  the  N2(A)  (or  He*)  inlet,  the  number 
density  of  N2(X,  v")  remains  essentially  unclumged.  This  is  because  vibrational  quenching  is 
a  relatively  slow  process  (k^  f  cm^  molecule'^  s*^ )  and  reagents  are  mixed  for  too 
short  of  time  (^2  ms)  to  effect  significant  quenching.  Any  diminution  in  the  first-positive 
(negative)  emission,  therefore,  results  from  quenching  of  N2(B)  (N2(B)),  or  pethaps  N2(A) 
(He*).  Injecting  the  reagent  at  the  upstream  end  of  the  flow  reactor  allows  adequate 
time  ( ~  10  to  30  ms)  for  the  N2(X,  v”)  to  be  quenched  unless  the  quenching  rate  coe^icients 
are  exceedii^ly  small.  Any  difference  in  dimirmtion  of  the  first-positive  (negative)  emission 
with  the  injector  in  the  upstream  position  as  compared  to  the  downstream  position  can  be 
attributed  to  quenching  of  N2(X,  v”)  by  the  added  reagent. 

Figure  4  illustrates  these  concepts.  It  slKiws  a  portion  of  the  spectrum  of  N2(B-A) 
emission,  excited  by  the  reaction  between  N2(A)  and  N2(^’  ^  absence  of  any  added 

quencher  and  in  the  presence  of  a  fixed  concentration  of  a  SiH4  at  a  small  mixing  time,  and 
one  21.7  ms  greater.  Adding  the  SiH4  at  the  short  time  drops  the  N2(B-A)  emission  about  a 
factor  of  1.6.  When  the  mixing  time  is  increased  by  21.7  ms,  the  first-positive  emission 
intensity  diminishes  almost  another  factor  of  1.6.  We  can  estimate  rate  coefficients  for 
N2(X,  v)  and  N2(B,  v)  quenching  by  SiH4  from  these  data  as  shown  below. 

In  the  absence  of  added  quencher,  the  intensity  of  the  N2  first-positive  emission  is 
given  by^* 


I^  =  k^2(B)l  =  kii(N2(X,v'')l[N2(A)l  .  (19) 


For  global  measurements,  k,,  tte  radiative  decay  rate  of  the  N2(B),^^  is  about  1.5  x  10^  s'^ 
TIk  experiments  use  a  helium  bath-gas  at  pressures  of  about  1  Torr.  This  reduces  quenching 
in  the  flow  tube  to  less  than  a  20%  effect.^^  As  a  result  we  neglect  bath-gas  quenching  in 
Eq.  (19). 

Then  the  quencher  is  introduced  into  the  flow  reactor  with  the  injector  in  the 
downstream  position,  that  is  short  interaction  times  where  N2(X,  v)  removal  is  negligible,  the 
first-positive  emission  intensity  is  given  by 
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Figure  4.  Variation  in  the  N2(B-A)  intensity  excited  in  the  N2(A)  +  N2(X,v)  reaction, 
as  a  function  of  reaction  time  and  SiH4  number  density.  Conditions  in  the 
order  of  decreasing  intensity  are:  a)  SiH4  -  0,  At  »  0  ms;  b)  SiH4  ~ 

2.3  X  10^^  molecules  cm*^.  At  =  0  ms;  c)  SiH4  =  2.3  x  10*^  molecules 
cm'^.  At  =  21.7  ms. 


I  =  k,lN2(B)l  = 


kn[N2(A)][N2(X,v'’)l 


(20) 


where  kg  is  the  rate  coefficient  for  N2(B)  quenching  by  SiH4.  N2(A)  number  densities  are 
also  measured  under  each  set  of  conditions  to  correct  any  quenching  of  the  N2(A).  Taking 
the  ratio  of  the  Hrst-positive  intensi^  in  the  absence  of  quencher  to  that  in  the  presence  of 
quencher  gives  the  classical  Stem-Volmer  formula^^  for  electronic  quenchii^,  viz.. 


VfNzAJo 

I/(N2A1 


=  1  +  kQ[Ca/k,  . 


(21) 


where  the  subscript  o  refers  to  conditions  in  the  absence  of  added  quencher.  The  rate 
coefficient  for  N2(B)  quenching  for  data  taken  with  the  injector  in  the  downstream  position 
is. 
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(22) 


With  the  injector  in  the  upstream  position,  the  N2(X,v'')  is  also  quenched  and  the 
first-positive  emission  intensity  is  describe  by 


I  *MN2(B)] 


kn[N2(X.v'')loe 

k. 


-MQlAt 


(23) 


where  k^  is  the  rate  coefficient  for  N2(X.  v")  quenchii^  and  At  is  the  time  the  quencher  and 
the  N2(X,  v")  are  mixed.  The  N2(X,  v")  quenching  rate  coefficient  is  determined  from  the 
slope  of  the  natural  log  of  the  ratio  to  Stem-Volmer  factors  at  long  to  short  decay  times 
plotted  against  [(^At,  i.e.. 


tn(Tjr^)  =  kJQlAt 


(24) 


where  the  subscripts  u  and  d  refer  to  measurements  with  the  injector  in  the  upstream  and 
downstream  positions,  respectively. 

5.2  Results  for  N2(X,  vSsS) 

Figures  5  to  7  are  representative  plots  of  the  data,  and  Table  1  summarizes  all  results 
from  the  studies  using  the  N2(A)  energy  transfer  diagnostic.  Since  rate  coefficients  for 
electronic  quenching  of  N2(B)  necessarily  are  determined  in  the  data  analysis,  they  are  also 
listed  in  Table  1.  The  error  bars  represent  one  standard  deviation  of  the  slopes  when  data 
are  plotted  according  to  Eq.  (24). 

The  vibrational  relaxation  rate  coefficients  determined  by  the  Eq.  (24)  analysis  should 
properly  be  taken  to  be  lower  limits  to  the  true  vibrational  relaxation  rate  coefficients.  Since 
the  1^3  mode  of  SiH4  is  resonant  to  within  kT  for  single  quantum  transfer  from  N2(X,v),  the 
mechanism  for  the  relaxation  most  likely  proceeds  one  step  at  a  time.  Thus  the  actual 
process  being  measured  will  be  the  net  difference  between  relaxation  of  all  vibrational  levels 
above  a  given  level  into  that  level  and  the  removal  of  vibrational  energy  from  the  given  level 
into  all  vibrational  levels  below  it.  Even  so,  the  lower-limit  vibrational  relaxation  rate 
coefficients  are  quite  lai^e,  being  about  twice  the  size  of  the  rate  coefficient  for  quenching 
N2(X,v)  by  N2O  and  an  order  of  magnitude  or  more  greater  than  the  rate  coefficients  for 
vibrational  relaxation  by  CO,  CO2,  and  N2,^^  species  that  are  generally  considered  to  be 
effective  at  removing  vibrational  energy  from  N2(X,v). 
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Figure  5.  Net  quenching  of  N2(X,  vs6)  by  silane 

The  rate  coefficients  for  electronic  quenching  of  N2(B)  seemed  unusually  large  and 
caused  some  concern  as  to  the  validity  of  our  measurements.  Since  the  primary 
measurement  is  of  the  ratio  of  two  intensities,  calibration  of  the  optical  system  will  not 
contribute  to  possible  errors  in  the  rate  coefficients.  The  other  two  measured  quantities, 
[SiH4]  and  At,  both  depend  primarily  on  gas  flow  rate  measurements.  Consequently  we  did 
an  extensive  set  of  calibrations  of  our  flow  meters.  These  calibrations  involved  measuring 
the  increase  of  pressure  in  a  calibrated  volume  as  a  function  of  time.  We  used  two  different 
volumes  that  we  calibrated  in  situ  by  comparison  with  a  flask  whose  volume  had  been 
determined  previously  by  measuring  the  difference  between  its  weight  with  the  flask  filled 
with  water  and  with  it  empty.  In  addition  we  used  two  different  pressure  transducers  for  the 
pressure  measurements.  The  results  of  these  measurements  showed  good  agreement  both 
with  previous  calibrations  and  with  manufacturer’s  specifications  for  a  number  different  gases 
and  flow  meters.  The  calibration  for  the  SiH4  flow  meter  was  tested  against  Ar,  N2,  CH4, 
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Figure  6.  Net  quenching  of  N2(X,  v^9)  by  SiH4 

N2P,  and  SFg  in  addition  to  SiH4.  With  the  exception  of  SiH4,  the  calibrations  were  within 
5%  of  manufacturer’s  specifications.  At  the  low  end  of  its  range,  the  SiH4  flow  meter  read 
about  15%  below  specifications.  At  the  top  end  of  the  flow  meter  range  the  discrepancy  was 
only  2%. 

We  also  considered  the  effects  of  N2(X,v)  quenching  during  the  short  transit  time 
from  the  SiH4  injector  to  the  observation  region  when  the  injector  was  in  the  downstream 
position.  This  correction  generally  amounted  to  no  more  than  10  to  20%.  Finally,  we  have 
used  this  same  technique  to  measure  rate  coefficients  for  N2(B)  quenching  by  a  variety  of 
other  gases  with  reasonable  agreement  between  our  results  and  those  of  other  groups  who 
used  different  measurement  techniques.  It  appears,  therefore,  that  silane  quenches  N2(B) 
with  remarkable  efficiency.  The  attractive  forces  between  the  two  molecules  are  sufficiently 
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Figure  7.  Net  quenching  by  N2(X,  v^ll)  by  SiH4 


strong  that  quenching  takes  place  at  a  rate  significantly  greater  than  the  hard-sphere  collision 
rate. 


5.3  Results  for  N2(X,  v^6) 

An  analysis  of  the  variations  in  the  number  densities  of  N2(X,v  <  6),  determined 
from  the  He*  Penning-ionization  diagnostic,  as  a  function  of  the  product  [SiHJAt, 
i.e.,  equivalent  to  the  analysis  related  to  Eq.  (24),  gives  plots  such  as  Figures  8  and  9  and 
the  set  of  rate  coefficients  given  in  the  second  column  of  Table  2.  As  in  the  results  in  Table 
1,  these  values  can  only  properly  be  considered  lower  limits. 

Because  the  Penning-ionization  diagnostic  provides  absolute  number  densities  of 
molecules  in  each  vibrational  level,  these  data  are  amenable  to  an  alternative  analysis  that 
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Table  1 .  Rate  Coefficients  for  N2(X,v)  and  N2(B^ng)  Quenching  by  SiH4 


N2(X,  v) 

N2(B,  V) 

Vibrational 

Levels 

Rate  Coefficient 
(10**^  cm^  molecule'^s’*) 

Vibrational 

Levels 

Rate  Coefficient 
(10'^  cm^  molecule'^ s' ^) 

^5 

4.1 

±0.5 

1 

2.2 

±0.5 

5.7 

±0.3 

2 

2.7 

±0.6 

5.3 

±0.4 

3 

2.2 

±0.5 

5.7 

±0.3 

4 

2.2 

±0.7 

5.9 

±0.4 

5 

2.8 

±0.8 

^9 

6.0 

±0.4 

6 

2.9 

±0.9 

^10 

5.9 

±0.5 

7 

2.9 

±1.2 

^11 

5.5 

±0.6 

8 

3.4 

±0.9 

^11 

6.4 

±0.5 

9 

3.4 

±0.8 

^12 

5.1 

±0.7 

4.3 

±1.3 

^13 

5.2 

±0.8 

11 

3.0 

±0.8 

213 

4.9 

±1.1 

12 

2.7 

±1.7 

Unweighted  avg 

5.5  ±  0.6 

Unweighted  avg 

2.9  ±0.6 

incorporates  the  effects  of  single  quantum  transfer  from  a  given  level  to  the  one  below  it. 
Although  the  form  of  our  data  set  makes  the  single  quantum  analysis  somewhat  crude,  the 
results  are  still  enlightening. 

The  differential  equation  describing  the  time  rate  of  change  of  molecules  in  level  v 
with  number  density  is 


dN^ 

“dT 


=  k^.iN^.i[SiH4]  -  k^N,[SiH4] 


(25) 


where  the  first  term  on  the  right-hand  side  of  the  equation  represents  the  rate  of  formation  of 
level  V  from  level  v+1  above,  and  the  second  term  represents  the  removal  rate  from  level  v 
into  level  v-1.  Rearranging  this  equation  and  integrating  it  gives 

j  dN,  =  N,  -  N,°  =  k,.,[SiH4j  I  N,.,dl  -  lc,(SiH4]  |  N.dt  .  O-d 

An  approximate  value  for  the  integral  expressions  on  the  right-hand  side  of  Eq.  (26) 
is 


1 


Nydt 


2 


At  . 


(27) 
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Figure  8.  Net  quenching  of  N2(X,  v=2)  by  SiH4 
Thus  we  can  rewrite  Eq.  (27)  as 

N,  -  N,"  -  k,.,  [SiH^]  At  -  k,  [SiH,]  At  . 


.  Summing  the  equation  from  level  v  up  to  level  v^,  the  highest  vibrational  level 
observed,  gives 
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Figure  9.  Net  quenching  of  N2(X,  v=4)  by  SiH4 
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Since  the  vibrational  po(Rilation  falls  off  fairly  rapidly  as  v  increases,  the  term  ii 
can  be  neglected  compared  to  the  term  in  v  provided  v  is  several  quanta  below  v^ 
Thus  after  rearrangement  we  are  left  with  a  simple  two  term  expression: 


Table  2.  Rate  Coefficients  for  Relaxation  of  N2(X,v)  Using  Penning-ionization  Diagnostic 


Vibrational 

Level 

Rate  Coefficient* 

Net 

Removal 

Single  Quantum 
Transfer  Analysis 

1 

0.84  ±0.05 

1.3  ±0.1 

2 

1.8  ±0.1 

2.6  ±0.2 

3 

2.8  ±0.1 

3.6  ±0.2 

4 

3.8  ±0.1 

6.0  ±1.0 

5 

4.8  ±0.2 

6.0  ±1.0 

6 

5.8  ±0.2 

4.5  ±1.0 

^Rate  Coefficients  at  300  K  in  units  of  10' 

*^cm^  molecule’*s‘* 

N  ° 


=  ky  [SiH4]  At 


(30) 


What  this  equation  says  is  that  the  number  of  molecules  lost  from  level  v  is  not  just  those 
determined  from  the  direct  change  in  the  population  of  level  v  itself,  but  actually  includes  the 
sum  of  the  population  lost  by  all  molecules  in  vibrational  levels  above  v. 

Figures  10  and  11  show  some  of  our  data  plotted  according  to  Eq.  (30).  The  slopes 
of  the  lines  in  these  plots  give  the  value  of  the  vibrational  relaxation  rate  coefficients.  These 
are  tabulated  in  the  third  column  of  Table  2.  The  analysis  via  single  quantum  transfer  gives 
somewhat  larger  values  for  the  rate  coefficients  for  relaxation  of  v  =  1  and  2,  but  approach 
the  values  derived  from  an  analysis  via  Eq.  (24)  for  higher  vibrational  levels. 

The  least  squares  fits  to  the  plots  related  to  the  Eq.  (30)  analysis  have  near  zero 
intercepts  for  the  lowest  vibrational  levels,  but  the  intercepts  become  increasii^ly  larger  for 
the  higher  vibrational  levels.  This  probably  results  from  neglecting  the  peculations  of 
vibrational  levels  above  6  in  the  analysis.  As  the  N2  vibrational  distribution  becomes 
increasingly  more  quenched,  the  populations  initially  in  higher  levels  will  caset^e  down  and 
contribute  significantly  to  the  populations  of  the  lower  levels.  That  is,  when  very  little 
overall  quenching  has  occurred,  vibrational  levels  7  and  above  will  contribute  negligibly  to 
the  populations  of  levels  such  as  3  or  4.  When  quenching  becomes  more  signifreant, 
however,  much  of  the  residual  population  in  levels  such  as  3  and  4  may  well  have  originated 
in  levels  7  and  above.  Since  the  Penning-ionization  diagnostic  cannot  provide  information 
about  the  populations  of  these  higher  levels,  our  analysis  necessarily  excludes  them  from 
consideration.  Were  it  possible  to  include  then,  the  ordinates  of  the  points  at  the  larger 
values  of  [SiH4]At  would  undoubtedly  be  higher.  This  effect  would  both  increase  the  slope. 
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Figure  10.  Single  quantum  analysis  of  quenching  of  N2(X,v  =  1)  by  SiH4 

thus  iiipr^ising  the  rate  coefficient,  and  also  would  lower  the  intercept,  indicating  that  a 
sufficiently  large  set  of  vibrational  levels  bad  been  included  in  the  analysis. 

The  rate  coefficient  for  relaxing  vibrational  level  1  can  also  be  determined  from  an 
analysis  of  the  increase  in  pqjulation  of  vibrational  level  0.  In  this  case  the  analysis  equation 


,No  -  No  . 

N.  *  N.° 


k,.,[SiH4]  A« 


Figure  11.  Single  quantum  transfer  analysis  of  the  quenching  of  N2(X,v  =  3)  by  SiH4 

This  equation  results  because  v  »  0  has  no  loss  processes,  and  because  in  single 
quantum  transfer,  only  those  molecules  removed  from  v  «  1  can  enter  v  «  0.  Figure  12 
shows  the  data  plot  relevant  to  Eq.  (31).  The  rate  coefficient  determined  from  these  data  is 
1.2  X  lO*^^  cm^ molecule**  s**,  in  excellent  agreement  with  the  value  determined  by 
considerii^  the  decay  of  v  »  1. 
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Figure  12.  Analysis  of  the  quenchii^  of  N2(X,  v=l)  by  SiH4  from  observations  of  the 
increase  in  [N2(X,  v=0)I  as  a  function  of  [SiHJ 
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6.  SiH4(»3)  EXCITATION  BY  NjCX,  v) 


We  used  the  apparatus  shown  in  Figure  2  to  investigate  the  excitation  of  SiH4(i'3)  by 
vibration-vibration  transfer  from  N2(X,  v)  to  SiH4.  Because  the  of  SiH4  is  resonant 

within  kT  with  single  quantum  transfer  from  all  N2(X,  v^  15),  excitation  seemed  the  most 
logical  choice  for  the  efficient  quenching  reaction.  The  1*3  mode  fundamental  frequency  is 
2190  cm'*  (4570  nm)  and  its  bj^  strength  is  625  cm'^  atm'*  at  STP.^  We  compared  our 
observations  with  observations  N2C)(i'3)  excitation  by  N2(X,  v)  when  N2O  was  added  to  the 
reactor.  The  N2O  mode  radiates  at  2224  cm'*  (4500  nm)  and  its  band  strength  is 
1541  cm'2  atm'*  at  STP.^* 

Figure  13  shows  the  fluorescence  spectrum  between  4  and  5  iim  when  SiH4  is  added 
to  a  flow  of  N2(X,  v).  The  strong  Q-branch  of  the  SiH4  mode  is  visible  at  4.58  nm  with 
associated  P-  and  R-branches  to  either  side.  The  spectral  resolution  is  0.04  /xm,  and  the 
SiH4  number  density  is  about  1.5  x  10*^  molecules  cm'^.  For  comparison.  Figure  14  shows 
the  fluorescence  spectrum  excited  when  N2O  is  added  to  N2(X,  v).  Here  the  spectral 
resolution  is  0.01  /im  and  the  N2O  number  density  is  about  5  x  10*^  molecules  cm'^.  In 
addition  to  the  N20(i'3)  fundamental  band  centered  at  4.50  /tm,  N2O  hot-band  emission  is 
evident  from  p^^l  at  4.56  /xm  and  i'3=3  at  4.62  /xm.  This  demonstrates,  that  in  the  case  of 
N2O  quenching  at  least,  multiquantum  vibrational  energy  transfer  is  efficient. 

Figure  15  shows  how  the  emission  from  SiH4(i'3)  at  4.56  /xm  varies  as  a  function  of 
SiH4  number  density  added  to  the  reactor.  At  low  number  densities,  the  emission  intensity 


Figure  13.  Spectrum  of  SiH4  {p-^)  emission  excited  by  N2(v) 
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Figure  14.  Spcctnim  of  emission  excited  by  NjCv) 


Figure  15.  Intensity  of  SiH4  (y^)  as  a  function  of  [8^4]  added  to  N2(X,v) 
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increases  almost  linearly  with  increases  in  SiH4  number  density.  For  higher  SiH4  additions, 
however,  the  emission  intensity  increase  is  less  than  linear.  Two  different  effects  are 
responsible  for  this  decrease  in  linearity.  The  higher  SiH4  number  densities  reduce  the 
effective  N2(X,v)  number  density  in  the  detector’s  field  of  view,  tlwreby  lowering  the  rate  of 
SiH4(i'3)  excitation.  In  addition,  higher  number  densities  of  added  SiH4  quench  the  SiH4 


emission.  We  can  test  these  suppositions  with  a  simple  kinetic  model. 

We  base  this  model  upon  the  following  set  of  reactions: 

SiH4  +  NjfX.v) 

-  SiH4(w)  -1-  N2(X.v-1) 

(32a) 

-*  other  products 

(32b) 

SiH4(w) 

-►  SiH4  +  hi/ 

(33) 

SiH4(w)  -»■  SiH4 

-*  2SiH4  . 

(34) 

SiH4(w)  -¥■  Ar,  N2 

-*  SiH4  +  Ar,  N2 

(35) 

Although  reaction  (35),  bath-gas  quenching,  could  not  be  deduced  from  the  data  in 
Figure  15,  the  necessity  for  its  inclusion  becomes  apparent  when  SiH^  excitation  is  studied  at 
several  different  pressures. 

With  the  above  reaction  scheme,  we  can  write  down  and  solve  two  differential 
equations  describing  the  temporal  variations  in  the  number  densities  of  N2(X,v)  and  SiH4(w). 
For  N2(X,v)  we  have 


d[N2(X,v)] 

dt 


-k32lSiH4][N2(X,v)] 


(36) 


which  can  be  solved  to  give 

[N2(X,v)]  =  [N2(X.v)1„  exp(-k32(SiH4]t)  .  (37) 


where  [N2(X,v)1q  is  the  N2(X,v)  number  density  in  the  absence  of  SiH4.  For  SiH4(w)  we 
have 


d[SiH4(w)] 

dt 


k32a[SiH4]IN2(X.V)l  -  (k33*k34[SiH4] 


^  kjsIAr,  N2])[SiH4(w)] 


(38) 
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Using  Eq.  (37)  for  [N2(X,  v)],  Eq.  (38)  can  be  solved  to  give 


[SiH4(w)]  = 


k32alN2(X,  (exp( -kjjlSiHJi)  -  exp( -(kjj  ♦k34lSiH4J  ♦  kjjlAr.  Nj])!)) 

kja-kjJSiHJ  -  kjjIAf.  NjJ  -  kjjtSiHJ 


Because  we  are  looking  along  the  axis  of  the  flow  reactor,  our  effective  field  of  view 
encompasses  emissions  from  volume  elements  at  all  axial  positions,  and  thereby  reaction 
times,  between  zero  and  the  bend  in  the  flow  reactor,  36  cm  distant  from  the  SiH4  injector. 
Thus  the  observed  emission  is  in  reality  the  integral  over  time  of  emission  from  each  of  these 
volume  elements.  .The  product  of  the  collection  efficiency  times  the  area  of  each  volume 
element  normal  to  the  flow  mbe  axis  is  constant  until  the  monochromator  field  of  view 
becomes  comparable  to  the  cross-sectional  area  of  the  flow  tube.  This  happens  in  our  reactor 
about  24  cm  downstream  from  the  injector.  Over  the  last  12  cm  of  the  field  of  view,  the 
product  of  collection  efficiency  times  volume  decreases  slowly  so  that  at  the  farthest  reaches 
of  the  flow  reactor,  the  effective  photon  emission  rates  are  attenuated  by  almost  a  factor  of  2 
compared  to  those  arising  from  emission  in  the  first  24  cm  of  the  reactor.  Rather  than 
account  for  this  effect  exactly,  we  found  it  convenient  to  correct  our  observations  by 
assuming  constant  collection  efficiency  over  a  column  shorter  than  actually  observed.  This 
allows  for  easier  integration  of  Eq.  (39)  over  time  and  has  a  negligible  effect  on  the  results 
provided  the  effective  length  of  the  coliunn  is  chosen  correctly. 

The  observed  emission  intensity  is  then  the  product  of  k33  times  the  integrated  number 
density  of  SiH4(w)  along  the  length  of  the  observation  volume: 


11f32ak33lN2(X.v)]JSiH4J 
k33*k34[SiH4]  >  kjjlAr,  Njl-  k32lSiH4J 


(l-exp(-k32lSiH4]to))  (l-exp('(k33+k34[SiH4)  +  k35lAr, 

k32lSiH4]  "  (k33*k34[SiH4l  ^  k35lAr.  I^2l) 


(40) 


where  t^  is  the  effective  time  along  the  length  of  the  observation  region  and  the  factor  t  is  an 
efficiency  factor  for  the  detection  system. 


We  subjected  the  data  in  Figure  15  to  a  least-squares  analysis  according  to  Eq.  (40) 
with  the  unknown  parameters  being  the  product  ik32j|k33[N2(X,v)]Q  and  k^^.  We  determined 
a  value  for  k32  from  our  earlier  work  to  be  between  1  and  5.5  x  10"^^  cmr  molecule'^  s*^ 
depending  upon  vibrational  level,  and  k33  can  be  calculated  from  the  measured  SiH4  1*3  band 
stren^.  Given  a  band  strength  of  625  cm'^  atm'^  we  calculate  a  radiative  decay  rate  of 
28  s'^  for  the  triply  degenerate  SiH4(i'3).  The  results  of  the  least-squares  fit,  for  an  assumed 
value  k32=2.5  x  10*^^  cm^  molecule"'  s"',  are  3.1  x  10""  for  the  product  {k32a*^33[N2(X,v)]Q 
and  4.5  x  10"'^  cm^  molecule"'  s"'  for  k34.  The  solid  line  in  Figure  16  compares  the  results 
of  the  fit  to  the  data. 
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Figure  16.  Intensity  of  SiH4(v3)  emission  as  a  function  of  [SiH4]  added  to  N2(X,v).  The 
solid  line  is  a  least-squares  fit  to  Eq.  (40)  with  k32  =  2.5  x  10'*^  cm^ 
molecule"^  s'^  and  k33  =  28  s 

To  compare  with  the  SiH4  excitation  data,  we  studied  the  excitation  of  N20(('3)  under 
identical  conditions.  This  affords  comparison  of  the  relative  excitation  efficiencies  between 
the  two  reactions.  A  set  of  equations  analagous  to  Eqs.  (32)  through  (40)  also  holds  for  the 
N2O  data  analysis.  When  necessary  we  shall  prime  the  k’s  for  the  rate  coefficients  for  the 
N2O  reactions  to  distinguish  them  from  those  for  the  SiH4  reactions. 

Data  taken  at  several  different  pressures  are  shown  in  Figure  17  for  SiH4  excitation 
and  Figure  18  for  N2O  excitation.  It  is  clear  from  the  way  the  initial  slopes  of  the  SiH4  data 
become  much  smaller  for  the  higher  pressure  data,  in  contrast  to  the  near  constancy  for  the 
N2O  data,  that  bath  gas  quenching  of  SiH4(i'3)  is  fairly  significant.  The  variations  in  the 
initial  slopes  can  result  fimm  the  effects  of  reaction  35  or  35'  or  from  changes  in 
[N2(X,  v)]p.  Because  at  each  pressure  the  value  of  (N2(X,  v)],,  will  be  the  same  for  both 
sets  of  data,  we  should  be  able  to  analyze  the  data  to  allow  estimation  of  k35.  In  addition, 
by  comparing  the  observed  intensities,  we  should  be  able  to  determine  a  ratio  of  ■ 

Values  for  k32,  k32',  and  k33,  k33'  can  be  estimated  from  the  quenching  data  and  band 
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Figure  17.  Excitation  of  SiH4(i>3)  by  N2(X,  v)  at  various  total  argon  pressures 
compared  to  a  kinetic  model 

strength  measurements,  respectively.  The  band  strength  data  give  values  of  k33  and  k33'  of 
28  and  214  s' ^  respectively,  while  our  quenching  data  indicate  a  value  for  k32  and  k32'  of 
about  2.5  X  lO'^^  cm^  molecule'^  s'^  airf  1  x  10'*^  cm^  molecule'^  s'*  respectively.  For  the 
two  parameters  we  can  determine  from  the  data  of  Figures  17  and  18,  the  actual  values  used 
for  k32  and  k32'  aren’t  particularly  important.  Variations  in  both  quantities  of  30  to  50% 
give  similar  results  for  k35  and  the  ratio  The  particular  values  chosen  for  k34, 

k34'  and  k35'  can  be  varied  almost  at  will  without  affecting  the  flnal  results  for  for  the 
quantities  k32a/1^2a'  results  indicate  the  ratio  ^  ^3S 

is  (6 ±2)  X  10'*^cm^  molecule'*  s'*. 

Given  that  the  ratio  k32/k32'  is  about  2.5,  the  result  above  indicates  that  N2(X,  v) 
energy  transfer  to  N2O  is  about  7  times  more  efficient  at  producing  infrared  radiation  at 
4500  nm  than  is  energy  transfer  to  SiH4.  Put  another  way,  if  the  branching  ratio, 
is  unity,  then  only  about  15%  of  the  time  will  an  N2(X,  v)  quantum  transferred  to  SiH4 
appear  as  radiation  at  4500  nm.  One  reason  for  this  discrepancy  is  undoubtedly  related  to 
the  fact  that  SiH4  has  an  infrared  inactive  vibrational  mode,  v^,  that  is  also  resonant  with 
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IN201(10'2  motecules  cm-3) 

Figure  18.  Excitation  of  N20(i'3)  by  NjfX,  v)  at  various  total  pressures  of  argon 
compared  to  a  kinetic  model 

N2(X,  v)  one  quantum  transfer.  Evidence  for  excitation  of  this  band  could  appear  in  the 
infrared  in  a  combination  band.  Although  we  searched  the  region  between  2000  and 
5000  nm,  we  saw  no  other  emissions  that  we  could  ascribe  to  SiH4  infrared  radiation.  Our 
system  is  not  particularly  sensitive,  however,  and  these  bands  are  likely  to  be  much  weaker 
than  the  band. 
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7.  ALTERNATIVE  DETERMINATION  OF  NjCX,  ▼)  QUENCHING 

RATE  COEFFICIENTS 


We  measured  rate  coefficients  for  N2(X,  v)  removal  by  SiH4  by  monitoring  the 
attenuation  in  N20(('3)  radiation  when  SiH4  was  added  to  a  flow  of  N2(X,  v)  in  a  SO-cm- 
long,  2.5-cm-diameter  flow  tube  that  was  attached  at  its  downstream  end  to  the  5-cm- 
diameter  flow  tube  used  to  monitor  the  infrared  radiation.  The  N2O  was  added  through  the 
loop  injector  in  the  S-cm-diameter  tube  directly  in  front  of  the  CaF2  window.  In  this  way, 
N2O  is  used  as  a  tracer  of  N2(X,  v)  and  changes  in  N20(i'3)  intensity  can  be  related  to 
changes  in  N2(X,  v)  number  density.  Unfortunately,  because  of  the  need  for  the  long  path 
length  to  obtain  adequate  signal  levels,  the  analysis  is  not  completely  straightforward. 

Proper  analysis  does,  however,  lead  to  reasonable  results. 

Three  effects  which  lead  to  a  reduction  in  N20(t'3)  intensity  when  SiH4  is  added  at  the 
upstream  end  of  the  2.5-cm-diameter  section: 

1 .  quenching  of  N2(X,  v)  in  the  2.5-cm  section; 

2.  quenching  of  N2(X,  v)  in  the  5-cm  section; 

3.  quenching  of  N2O  fluorescence  by  SiH4. 

The  expression  for  [N20(i»3)]  in  the  5-cm-diameter  section  as  a  function  of  distance  (time) 
along  the  length  of  the  field  of  view  is  similar  to  that  given  above  in  Eq.  (39).  One  major 
difference  is  that  the  factor  (N2(X,  v)]o  in  Eq.  (39)  becomes  [N2(X,  v)](,  exp(-k32{SiH4]At,), 
where  Ati  is  the  transit  time  in  the  2.5-cm-diameter  tube  between  the  SiH4  injector  and  the 
entrance  of  the  5-cm-diameter  tube.  Since  our  observations  are  actually  proportional  to  the 
integrated  fluorescence  intensity  along  the  length  of  the  5-cm-diameter  section,  we  have  to 
integrate  the  product  of  k33'[N20(i'3)]  over  the  time  the  reactants  are  in  the  field  of  view, 

At2.  The  resultant  expression  is 


k32k„ [N2O]  [N2(X,  v)1„C  -'‘nCSWJAt, 


1  .  g-(knIN20)»k32[SiHJ)At2 


k33+k32lSiH4]  -  (k„[N201  +  k32lSiH4]) 


k;2[N20]  *  k32[SiH4] 


k„  *  k34[SiH4] 


(41) 


For  the  experiments  we  ran,  this  expression  can  be  simplified  because  the  product 
k33'At2  is  sufficiently  large  that  the  second  exponential  term  can  be  neglected.  After  this 
simplification  and  some  algebraic  manipulation,  the  expression  for  the  observed  intensity 
becomes 
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k^lNzO]  (N2(X.  v)1„  C  ■‘‘“fSiHJAt, 

(-kjjfNjOI  * 

1  _  ® 

(42) 

kjjlNjO]  *  k32[SiH4] 

1  + 

.  *^33 

hi 

1  ♦ 

^33 

We  fit  several  sets  of  decay  data  to  the  natural  log  of  Eq.  (42),  using  explicit  values 
for  the  rate  coefficient  k32'.  Figure  19  compares  several  sets  of  data  with  the  best  fit  curves 
described  by  Eq.  (42).  The  results  of  the  fits  wve  k32  =  (2.5  ±  0.5)  x  10'^^  cm^ 
molecule'^  s‘*,  and  k34'  =  (1.1  ±  0.4)  x  10*"  cm^  molecule**  s**.  These  values  were 
insensitive  to  50%  variations  in  k32'.  The  value  for  derived  from  this  analysis  is 
consistent  with  values  derived  by  monitoring  decays  in  N2(X,  v).  Because  the  fluorescence 
of  N2O  can  be  excited  by  quenching  a  number  of  different  vibrational  levels  of  nitrogen,  we 


Figure  19.  Quenching  of  N2(X,  v)  by  SiH4  using  J42*^(*'3)  fluorescence  as  a 
monitor  of  [N2(X,  v)] 
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would  expect  this  analysis  would  give  a  value  for  somewhat  larger  than  that  for 
quenching  vs i,  but  smaller  than  we  had  observed  for  quenching  v^S  or  6.  This  is  exactly 
what  is  observed. 
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8.  SILICON  NITRIDE  COATING  EXFERlMEr^ 


We  investigated  briefly  the  coating  of  sapphire  and  single-crystal  silicon  wafers  with 
silicon  nitride.  We  used  the  technique  of  infrar^  absorption  spectroscopy  to  analyze  the 
coatings  we  produced.  In  addition  to  coating  the  substrates,  the  experiments  generated 
copious  quantities  of  white  powder.  Infrared  absorption  spectra  of  KBr  pellets  with  small 
quantities  of  this  powder  mixed  in  showed  the  powder  to  be  silicon  nitride  also.  Figures  20 
and  21  show  typical  spectra  of  the  film  on  a  silicon  substrate  and  of  the  powder, 
respectively,  llie  strong  band  at  around  900  cm'^  shows  that  the  samples  are  silicon  nitride. 
The  presence  of  bands  around  3300  cm‘^  and  1400  to  1050  cm'^  indicate  the  presence  of 
N-H  and  N-H2  bonds  in  the  samples,  while  the  small  band  at  2150  cm'^  shows  that  some  SiH 
bonds  are  also  present. 

The  structure  in  the  spectrum  in  Figure  20  between  1800  and  1400  cm'^  is  not  a 
property  of  the  coating,  but  rather  a  peculiarity  of  the  single-crystal  silicon  wafers  used  in 
the  coating  process.  The  spectrum  was  taken  with  an  uncoated  wafer  in  the  reference  beam 
of  the  spectrophotometer.  The  spectrum  of  an  uncoated  wafer  between  2000  and  1200  cm'^ 
is  also  shown  in  Figure  20.  The  uncoated  wafer  also  shows  strucmre  in  the  1800-1400  cm'^ 
region.  Apparently  the  small  strucmres  in  the  substrate  spectra  are  amplified  when  a  wafer 
is  used  in  the  reference  beam.  The  referencing  procedure  is  important,  however,  since 
silicon  has  several  small  absorption  bands  in  the  spectral  region  covered.  These  bands  are 
removed  from  the  spectral  signamre  of  the  coating  by  using  an  uncoated  wafer  in  the 
reference  beam. 
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Figure  20.  Infrared  absorption  spectrum  of  Si3N4  on  single-crystal  silicon 
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Wavelength  in  bun) 
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Microscopic  examination  showed  the  central  portions  of  the  coatings  to  be  relatively 
smooth  and  uniform.  The  coating  could  be  scndched  rather  easily,  however.  We  suspect 
this  results  from  the  incorporation  of  hydrogen  into  the  coating.  We  think  that  reducing  the 
flow  rate  of  silane  relative  to  the  active  nitrogen  would  probably  lead  to  reduced  hydrogen 
incorporation.  Unfortunately,  we  were  thwarted  in  our  attempts  to  test  this  hypothesis. 

Our  initial  apparatus  was  limited  in  terms  of  the  minimum  SiH4  flow  rates  we  could 
achieve.  The  relatively  high  SiH4  flow  rates  prevented  conq)lete  decomposition  of  the  silane, 
which  probably  explains  the  presence  of  so  much  hydrogen  in  the  films  and  powder.  In 
order  to  have  smaller  flows  of  SiH4,  we  made  a  5%  mixture  of  SiH4  in  argon  and  stored  it 
in  a  5  1  bulb  on  our  vacuum  line.  Several  coating  experiments  with  this  mixture  were 
disappointii^.  We  apparently  developed  an  air  leak  in  the  gas  handling  system  because  our 
coatings  with  the  lower  Sill4  proved  to  be  silicon  oxide  rather  than  silicon  nitride. 

The  silicon  nitride  powder  was  generally  very  fine  and  had  a  tendency  towards  easy 
agglomeration.  It  was  also  very  light  and  was  easily  dispearsed  by  the  merest  of  air 
disturbances.  In  addition,  it  was  strongly  affected  by  sinall  static  charges  on  the  surfaces  of 
various  laboratory  components  such  as  sample  jars  and  microscope  slides.  Microscopic 
observation  showed  that  the  individual  particles  were  on  the  ordn  of  a  micron  or  less  in 
diameter.  In  a  few  instances  the  powdn  adhered  more  tightly  to  the  walls  of  the  flow 
reactor  and  was  removed  essentially  as  small  flakes.  The  KBr  pellet  spectrum  of  the  flakes 
was  essentially  the  same  as  the  fine  powder.  Qualitatively,  it  aj^ieared  that  the  quantity  of 
powder  formed  correlated  with  the  N-atom  number  densities  in  the  active  nitrogen.  Although 
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we  did  implement  a  quantitative  monitor  of  N-atom  mimber  densities,  it  was  too  late  in  the 
program  to  investigate  this  correlation  in  a  quantitative  fashion. 

Although  our  initial  kinetic  measurements  on  the  quenching  of  N2(X,  v)  by  SiH4  were 
quite  extensive,  we  saw  little  evidence  that  silane  was  decomposed  to  a  significant  extent  in 
these  experiments,  with  one  slight  exception.  The  walls  of  the  flow  reactor  used  in  the 
kinetic  studies  did  develq)  a  white  film  over  time  just  downstream  from  the  N2(A)  injector. 
Although  we  dkln’t  analyze  this  film  it  is  undoubt^y  a  coating  of  silicon  nitride  powder. 

The  reactor  walls  were  still  relatively  clear  upstream  from  the  injector.  It  would  seem, 
therefore,  that  the  N2(X.  v)  itself  is  not  sufficient  to  decompose  the  silane  sufficiently  to 
form  silicon  nitride.  The  presence  of  the  metastables,  on  the  other  hand,  did  allow  silicon 
nitride  fomuition  to  occur.  It  is  not  clear  whether  the  reaction  between  the  metastables  and 
silane  dissociated  the  silane  into  fragments  that  could  then  react  with  N2(X,  v)  to  form 
silicon  nitride,  or  if  the  silicon  nitride  formation  resulted  solely  from  tlw  interaction  of  the 
metastabfes  with  the  silane  and  its  subsequent  decomposition  fragments.  This  is  an  area  that 
merits  (urther  investigation. 


9.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  STUDY 


We  showed  that  silane  is  extremely  efficient  at  quenching  vibrational  energy  of 
N2(X,  v).  Only  about  15%  of  the  energy  transfered  results  in  excitation  of  SiH4(i'3).  The 
rest  presumably  excites  the  resonant,  but  non-infrared  active  mode  and  other  modes 
outside  the  spectral  bandpass  of  our  detection  system. 

The  success  of  the  single-quantum  transfer  analysis  indicates  that  single  quantum 
vibration-vibration  transfer  is  the  most  likely  ptuhway  for  the  energy  flow  in  the  system. 

Our  kinetic  observations  could  not  rule  out,  however,  the  possiblilty  of  multiquantum  transfer 
resulting  in  hot  band  excitation  of  SiH4  or  even  in  direct  SiH4  dissociation.  We  had  hoped  to 
be  able  to  monitor  H-atom  production  in  this  system,  but  a  number  of  equipment  failures 
prevented  our  being  able  to  follow  these  investigations. 

Our  limited  coating  studies  showed  that  adding  silane  to  an  active  nitrogen  afterglow 
resulted  in  facile  production  of  both  silicon  nitride  powder  as  well  as  in  coatings  on  sapphire 
and  silicon  substrates.  The  quantity  of  the  powder  produced  seemed  to  be  proportional  to  the 
number  density  of  atoms  in  the  afterglow,  but  we  were  unable  to  make  these  observations 
quantitative.  The  lack  of  significant  silicon  nitride  production  in  the  apparatus  used  for  the 
kinetic  measurements,  except  in  the  immediate  vicinity  of  the  metastable  injector,  shows  that 
vibrationally  excited  nitrogen  by  itself  does  iMt  appear  to  be  sufficient  for  silicon  nitride 
formation.  Whether  or  not  it  enhances  silicon  nitride  production  in  the  presence  of  atoms  or 
metastables  is  still  not  clear.  Since  the  reaction  of  atomic  nitrogen  with  silane  is  endoergic, 
it  seems  likely  that  N2(X,  v)  is  necessary  to  promote  efficient  silicon  nitride  synthesis.  This 
may  not  be  the  case;  however,  it  could  be  adequate  just  to  have  N  atoms  and  electronically 
excited  molecular  metastables  present.  The  primary  product  of  N-atom  recombination  is 
N2(B  ^Ilg),  and  our  kinetic  studies  show  that  silane  quenches  this  species  with  extraordinary 
efficiency. 

Future  investigations  should  focus  more  directly  on  the  coating  process  itself.  Using 
techniques  we  have  developed  and  implemented  in  this  program,  measuring  coating 
efficiencies  under  conditions  where  the  the  N-atom  and  N2(X,  v)  number  densities  are 
carefully  monitored  and  controlled  should  be  straightforward.  If  N2(X,  v)  appears  not  to  be 
particularly  important  in  the  silicon  nitride  formation  processes,  then  the  relative  importance 
of  electronically  excited  nitrogen  molecules  in  the  silane  decomposition  processes  should  be 
investigated.  This  can  be  done  in  two  ways.  Varying  the  reactor  pressure  under  conditions 
of  constant  N-atom  number  density  allows  the  production  rates  of  electronically  excited 
nitrogen  to  be  varied  in  a  controlled  fashion.  In  addition,  one  could  make  metastables  in  the 
absence  of  atomic  nitrogen  and  monitor  changes  in  silicon  nitride  formation  when  these 
metastables  are  added  to  a  flow  of  N  atoms. 
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